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Abstract. We describe an operation on dimer configurations on the hexagon lattice, called 
"squishing" , and use this operation to explain some of the properties of dimer generating 
functions. 

1. Introduction 

In this paper, we will describe and use a novel technique called "squishing", which one 
applies to the dimer model on the regular honeycomb lattice. We developed this technique 
as an attempt to verify a conjectured generating function which arises in algebraic geometry 
(specifically, in the Donaldson-Thomas theory of the orbifold C 3 /Z2 x Z 2 [lj). Our attempt 
was only partially successful, and we were later able to compute the generating function by 
other means. However, the technique is interesting in of itself, being one of relatively few 
dimer model techniques which exploits the self-similarity of the lattice at different scales. 

We begin by describing the original motivation for this work. 

Definition 1. A 3D Young diagram (or 3D diagram,) n is a subset of (Z> ) 3 such that if 
(x, y, z) G n, then (x', y', z') 6 n whenever x' < x, y' < y, and z' < z. 

3D Young diagrams are called boxed plane partitions or 3D partitions elsewhere in the 
literature. We refer to the points in n as boxes - the point (i,j, k) corresponds to the unit 
cube with vertices {{i ± ^, j ± ^, k ± |)}. 

We will be discussing the following generating functions: 

Definition 2. A weighting of (Z> ) 3 is a map 

w : (Z> ) 3 -> {p,q,r, s}, 

where p, q, r, s are formal indeterminates. w(p) is called weight of the lattice point p; the 
weight of a 3D diagram is the product of the weights of the lattice points at the centers of all 
of its boxes. The ^-partition function is then defined to be the formal sum 

z w = 2J w ( n )- 

7r 3D diagram 

The weightings we will be concerned with are the monochromatic weighting tfm, 
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and the Z 2 x Z 2 weighting wz 2 xz 2 , 
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along with various specializations of these weightings; we shall denote their partition func- 
tions Z{xy(p) and Zz 2X z 2 (p, q, r, s), respectively. 
It is a classical result jl] that 

(1) J> {1} (7r) = M(l,p) 

7T 

where we define 

n=l v 7 

Equation ([I]) also arises in algebraic geometry, essentially because 3D Young diagrams are 
the same as monomial ideals / C C 3 [a;, y, z] (simply read off the exponents of the elements 
of the coordinate ring C 3 //; these are the boxes of it). As a result, ([I]) is, in a certain sense, 
an invariant of C 3 ; specifically, it is the Donaldson-Thomas partition function for the space 
C 3 , up to a sign on p. [S]. 

It is possible [Tj to develop Donaldson-Thomas theory for orbifolds of C 3 under the actions 
of certain finite abelian groups. It turns out that for the group Z 2 xZ 2 , the partition function 
is given by 



M(l,Q)*M(qr,Q)M(q S ,Q)M(rs,Q) 



M(-q, Q)M(-r, Q)M(-s, Q)M(-qrs, Q) 



where Q = pqrs and M(a, z) = M(a, z)M(a~ 1 1 z). 

The curious identity §2§ was proven in pQ using vertex operators; it was conjectured earlier 
by Jim Bryan (based on the behaviour of related Donaldson-Thomas partition functions) and 
independently by Rick Kenyon (in an equivalent form, based on empirical computer work). 

It is possible to check some of the properties of [2] in an elementary manner. For example, 
it should be the case that specializing p = q = r = s should give 

Z Z2XZ . 2 (p,p,p,p) = Z {1} , 

and indeed ^ does satisfy this relation. Another striking relation suggested by ^ is: 

(3) Z ZaxZ3 (p, -1,-1,-1) = M(1,Q) 2 ; 

however, the combinatorial reason for this is far less clear. This paper demonstrates ^ in 
an elementary manner, without relying on the theorems of [TJ. 



2. Matchings, and Weightings 

Our attack on ^ immediately requires us to to encode the "surface" of a 3D Young 
diagram with dimers on the hexagon lattice. Let us fix some terminology. 

Definition 3. Let G = (V, E) be a graph. A matching or 1-factor of G is a subgraph 
M = (y,E') such that the degree of every vertex v € V is 1. Equivalently, a matching of G 
is a partition of the vertices of G into a disjoint union of dimers, or pairs of vertices joined 
by a single edge of G. 
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Figure 1. A matching of #5,4,3. 
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Figure 2. A 3D Young diagram viewed as a matching on H : 



3,3,3- 





For a general survey of results on the dimer model, see [2]. 

We will be considering matchings on the semiregular hexagonal mesh of side lengths 
a,b,c,a,b,c, denoted H a ^, c (see Figure [T] for a definition-by-picture). Matchings on this 
graph are in bijection with 3D Young diagram which are contained within an a x b x c 
box. To see why, imagine viewing a 3D Young diagram from a distance(i.e., under isometric 
projection). The faces of the 3D diagram are then rhombi, each of which is composed of 
two equilateral triangles (see Figure [2]). The centers of these triangles fall at the vertices of 
H a ,b,c, so we get a matching by replacing each rhombus with the corresponding edge. 

A weighting of a graph assigns a monomial to each of the graph's edges. Our first example, 
called a monochromatic weighting, is the one depicted in figure [3j which assigns a weight 
of one to all non-horizontal edges, and weight p l to horizontal edges which have i other 
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Figure 3. A monochromatic weighting on #4,4,4 




horizontal edges directly below them. We will adopt the convention that if an edge has no 
weight written beside it, then that edge has a weight of one. 

If our graph has a weighting, then the weight of a matching is the product of all of the 
weights of the edges of the graph. 

Now, we have defined two monochromatic weightings - one for 3D Young diagrams and 
one for graphs. One might ask whether they are "the same" : is the weight of a 3D diagram 
which fits inside an i x j x k box the same as the weight of the associated matching of Hi j^l 
Strictly speaking, the answer is no, because the weight of the empty 3D Young diagram 
should be zero, but the weight of the associated matching (see Figure [5]) is nonzero. 

However, given a matching M, we can define its normalized weight as the weight of M 
divided by the weight of the empty 3D Young diagram. Now it is easy to see that the 
normalized weight of M is equal to the weight of the associated 3D diagram ir. The proof is 
by induction on the number n of boxes in ir, the case n = being trivial. One only needs to 
check that the operation of adding a box to it has the effect of increasing both weights by p, 
which is easy to do. 

There are many other weightings on the hexagon meshes whose normalized versions are 
equivalent to the monochromatic weight. For example, one could rotate the weighted mesh 
by 120 degrees. Indeed, for our purposes, the following weighting is superior (see Figure [1]): 
we replace p with t, and superimpose three copies of the old monochromatic weighting: one 
normal, one rotated 120 degrees, and one weighted 240 degrees. This weighting assigns each 
box the weight t 3 , so substituting p = t 3 gives a new monochromatic weighting. 

Definition 4. The weighting described above is called w p . 

It is cumbersome to draw diagrams with p 1 ^ 3 as an edge weight, so throughout the paper, 
we shall use the convention that p = t 3 when it is convenient. 

We would next like to define a weighting of # a ,6,c whose normalized version is equivalent 
to the Z 2 x Z 2 weighting. There are several ways of doing this, but for our purposes, the best 
way is to first define three weightings: the qt-, rt-, and st-weightings (see Figure [6]). The 
gt-weighting is equivalent to the Z 2 x Z 2 weighting under the specialization p i— ► t; r, s \— > 1, 
and similarly for the other two weightings. 

Having done this, we construct the Z 2 x Z 2 weighting by assigning each edge in H a ^ c the 
product of its qt,rt, and st-weights. This weights boxes colored q, r, s correctly, but each 
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Figure 4. A "better" monochromatic weighting, w p , where p = t 3 




Figure 5. The empty 3D Young diagram and its associated w p - weighted matching 



/ V 
/ V 

/ '/ > X V N, V X. 

; / / x v 

/ i x v 
/ \ f yv y / v 




1 t t 1 

i £ i 

i i 
i 



box in the p position gets the weight t 3 . So specializing t i— ► p 1 / 3 gives the Z2 x Z2 weighting 
(see Figure [7|). Observe that we have assigned weight 1 to all of the grey edges. 



Definition 5. We call the weighting of Figure^ 



w 



p,q,r,s ■ 



3. Overlaying pairs of matchings 

We were introduced to the ideas in this section by Kuo's beautiful paper on graphical 
condensation [3]. In the following, G will always be a bipartite graph. 

Suppose that we have two matchings M\, M 2 of G. If we overlay these two matchings on 
the vertex set V of G, we have a multigraph N in which each vertex has degree two, called 
a 2-f actor of G. This terminology is slightly nonstandard in graph theory - elsewhere in the 
literature, a 2-factor is usually a collection of closed loops and isolated edges (not doubled). 
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Figure 6. The qt-, rt-, and st-weightings on #4,4,4 . 




If the edge e occurs in both Mi and M 2 , then e occurs as a doubled edge in N. In this case, 
since the degree of both endpoints of e is two, e is a connected component in N. Conversely, 
all doubled edges in N must occur in both Mi and M 2 . If we disregard the doubled edges, 
the rest of N decomposes into a collection of disjoint closed paths. 

Conversely, one can split a 2-factor into two one-factors: 

Lemma 6. A 2-factor N may be partitioned into an ordered pair of matchings (Mi, M 2 ) in 
precisely 2*^ closed paths in N > distinct ways. 

Proof. Suppose we have a 2-factor iV of a bipartite graph G. We may obtain two matchings 
Mi and M 2 of G as follows: If e is doubled in N, then place e into both Mi and M 2 . If P is 
a closed path in N, select one of the edges in P and place it into Mi. Place the next edge 
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Figure 7. The Z 2 x Z 2 weighting. 
1m* 




Figure 8. Overlaying two matchings on H 3 S 3 . 




in the path into M 2 , and so forth. Since G is bipartite, the path P is of even length, so each 
vertex in P has degree 1 in both Mi and M 2 . 

There are 2 ways to divide P between M 1 and M 2 , so there are 2 # ^ closed paths 111 N ^ pairs of 
matchings Mi, M 2 which correspond to N. □ 
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Figure 9. The hexagonal mesh -£^4,4,4 being squished. 




4. Squishing 

Consider the hexagonal mesh with even side lengths, i?2a,26,2c- The leftmost diagram in 
Figure [9] shows a picture of i/4,4,4, with some of the edges colored grey. The grey edges come 
in sets of three, all of which are incident to one central vertex. I shall call these three edges 
a propeller. 

The other two diagrams in Figure [9] show what happens when the length of the edges 
in each propeller is decreased, while the other edges remain long. We call this procedure 
squishing. Of course, the length that we choose to draw the edges in the graph has no bearing 
of the structure of the graph, but when the propellers are quite small, the graph of ^4,4,4 
"looks like" the graph of ^2,2,2 with each edge doubled. 

We will denote this squishing operation by the symbol ip : 

^ : H 2a ,2b,2c \ {propellers} — ► H aM 

where ip sends each edge to its image after squishing. It is ofted useful to speak of squishing 
a set E of edges of H 2a ,2b,2c, and we shall also denote this operation by ip: 

i>{E) = |J ^(e)- 

e6_B\{propellers} 

Sometimes, given E 1 C H a ^ c , we will need to look for sets of edges E for which ip(E) = E'. 
Naturally, there are many such E, since if) ignores the propellers and is two-to-one on all 
other edges. However, for a given E' C H aj b, c , there is a "most relevant" preimage of E' 
under ip, which we shall call (f(E') C H 2a ,2b,2c, defined as follows: 

E := |J ^(e), 
<p{E') := Eq U {all propellers adjacent to E } 



In other words, ip(E') contains all edges which squish to edges of E', as well as all propellers 
incident to those edges (See Figure 10). It is clear that ip o <p(E') = E'. 
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FIGURE 10. The "unsquishing" map (p 





Figure 11. A matching on f/444 being squished to a 2-factor on -££2,2,2- 
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5. Three lemmas on weightings of squished graphs 
When one draws a matching on the graph of -£^4,4,4 before squishing, we get what looks 



very much like a 2-factor of -££2,2,2 (See Figure 11 ). We will quantify precisely how this occurs 



and use the effect to prove several facts about certain specializations of the Z2 x Z2 partition 
function. 

Let us examine one of the propellers more closely, and determine what possible configu- 



rations of its edges can appear in a perfect matching M, up to symmetry, (see Figure 12). 
We label the central vertex D, and the other vertices A, B, C. First, observe that preciesly 
one of the three "short" edges must be included in M, to give the central vertex degree one. 
Suppose that it is the short horizontal edge CD. In order for vertices A and B to have 
degree 1 as well, they must each have one incident long edge. Up to symmetry, there are 
three ways for this to occur: 

(1) All three edges are horizontal; 

(2) The edge incident to A is horizontal; the edge incident to B is diagonal; or 

(3) The edges incident to both A and B are diagonal. 

If we now identify the vertices A, B, C, D and ignore the short edges, we can see that case 
1 gives rise to a doubled edge, and cases 2 and 3 give rise to a path through the vertex. 
Therefore, collapsing all propellers in H 2a ,2b,2c to points transforms a matching on H 2a ,2b,2c 
into a 2-factor on H a ^ c . In other words, the squishing operation i/j induces a map, which we 
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Figure 12. All possible configurations of edges around a propeller in a perfect matching. 




Figure 13. The Z 2 x Z 2 weighting, under the specialization £, q, r, s i— ► — 1 
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will call \P: 

^ : {1-factors of H 2a ,2b,2c} — ► {2-factors of H a ^ c }. 

It is easy to see that \I> is surjective but not injective. For example, #2,2,2 has 20 perfect 
matchings, whereas #1,1,1 has only 3 2-factors. 

It is often necessary to find the set of all matchings on ip(E'). Strictly speaking, this set 

is (ty\ip(E>)) {E')-, but for simplicity of notation we shall just write this as ^I~ 1 (E'). 

We will need to use the Z 2 x Z 2 weighting under the specialization t,q,r,s 1— > —1, which 
we shall call u>_i for short (see Figure 13). 

Lemma 7. Let A be a 2- factor of H abc . Then 

^ ^ t/7 ( ^ab+bc+ca _ ^^{closed paths in A} 
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Figure 14. The unique configuration of edges which squishes to a doubled edge. 
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Figure 15. Lifting a closed path £j to <f(£j), and completing it to a matching. 

-l -l 




Proof. Let us suppose that A decomposes as the disjoint union of doubled edges ei, . . . , e n 
and closed loops £±, . . . ,£ m . Because the unsquishings <f(ei) and ip(ftj) are all pairwise non- 
adjacent in H2a,2b,2c, it is clear that 



(4) £ «;_ 1 0x)=n £ n E 

*=1 \A» matching on ip(ej) / 3=1 \/^ matching on 

Let us first consider 1 (e^) , the preimage of a doubled edge. There is only one possible 



matching on 0(e$) (see Figure 14); it has weight -1. 

Next, let us consider \l/~ 1 (£j), the preimage of a closed loop. Every term in ty~ l (£j) is a 
matching on which maps to £j under ip. We shall "lift" £j up to <f(£j), so that at each 
step we choose an edge e' G ip(£j) which maps to e under ■0. If we do this in such a way 
that no two of the e' are adjacent, then there is a unique way to complete the lifted walk to 
a matching on ip(£j), by choosing exactly one edge of each propeller to be in the matching 



(see Figure 15). Note that the u>_i-weight of the lifting (second image) is the same as the 
weight of the matching (third image) since all of the propellers have weight 1. 

Suppose that we walk along the lifted path in ip(£j), in the counterclockwise (positive) 
direction, and we step on the edges e[, e' 2 , ■ • • , e' k , e[, in order. We keep track of the product 
of the weights of the edges we have stepped on. 

At any point in our walk, there are essentially only four states that we can be in. We can 
be standing on an edge weighted 1 or -1, and we can be facing one of two types of propeller: 
~ or ~ . Let us label the states 1 through 4, as follows: 

(1) Standing on a weight 1 edge, facing a ~ propeller 

(2) Standing on a weight -1 edge, facing a / propeller 

(3) Standing on a weight 1 edge, facing a ~ propeller 

(4) Standing on a weight -1 edge, facing a \ propeller 

To get from one edge to another in the path, we must turn either left or right at the 
propeller we are facing, so we shall think of our walk as a sequence of Rs and Ls, read 
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right-to-left. For example, the walk in Figure 15, starting at the upper leftmost horizontal 
edge and proceeding counterclockwise, is LRRLLLLRLLRLLL. 

Our trick is to interpret R and L as the following state-transition matrices: 
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These matrices encode which states can be reached from which, and keep track of the 
total weight of all paths leading from one state to another. For example, at the upper-left 
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most horizontal edge in Figure [15j we are in state 4, and we are about to turn to the left. 
We could go into state 1, picking up a weight of 1 (as in the diagram); alternately, we could 
go into state 2, picking up a weight of -1. This is modeled by our state-transition matrices, 
since 
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By the same token, the coefficients of the matrix LR are the sums of the weights of all paths 
which start in state i , turn first right and then left, and end in state j. If we want to know 



the sum of the weights of the lifts of the closed loop in Figure 15, we should compute the 
matrix product LRRLLLLRLLRLLL. We must end in the same state as we start since our 
path is closed; furthermore, our starting state must be either 3 or 4, so the total number of 
paths is the sum of the (3,3) and (4,4) entries of this matrix product. 

Now, here is the punchline: It is easy to check that L = i? -1 . In any closed loop on the 
hexagonal lattice, there must be 6 more left turns than right turns, so the combined weights 
of all lifts of any closed loop is the sum of the (3,3) and (4,4) entries of L 6 . Now, it is 
also easy to check that L 6 = —I, so this combined weight is -2, regardless of the shape or 
orientation of the loop £j. 

Going back to Equation [4j and recalling that there are i doubled edges and j closed loops 
in the 2-factor A, we can now state that 



W -M 



'_ 2 r 



'-!)*+' ■ 2 j 



Observe that i+j is the total number of connected components in A, so the proof is complete 
modulo the following lemma. □ 

Lemma 8. Let X be a 2-factor of H a ^ c . LetC(X) denote the number of connected components 
of A. Then C(A) = ab + be + ca (mod 2). 

Proof. Decompose A into two matchings, Ai and A2. Let us consider Ai for the moment. 

When one removes (or adds) a cube to the 3D diagram corresponding to Ai, one obtains 
a new matching. Graph-theoretically, this is equivalent to performing the change depicted 



in Figure 16 Let us call such an operation a r -transformation. Since it is possible to change 

any 3D diagram into any other by adding and removing boxes one at a time, repeated 

r- transfer mat ions will allow us to change any matching into any other matching. 
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FIGURE 16. A r-transformaion 
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FIGURE 17. Changing one 2-factor into another with r-transformations 



Let us alter the matchings Ai and A2 in this fashion, so that after finitely many steps, we 
have obtained two copies of the empty 3D Young diagram. At each stage, superimposing 
the two matchings gives a 2-factor, so we may think of r-transformations as acting on the 



set of 2-factors as well (see Figure 17) 



Next, we show that r-transformations leave C(A) invariant mod 2. Suppose we are apply- 
ing the r transformation to Ai, and we are affecting the edges around a particular hexagon 
H . Remove all edges of H from Ai to obtain a new subgraph A° of H at b iC . Now superimpose 
A^ and A2. The resulting graph has degree 2 everywhere except at vertices of H, all of which 
have degree 1. As a result, every vertex of H is connected to precisely one other vertex of H 
by some path (possibly of length one). Up to symmetry, there are two possible ways in which 



this may occur (see Figure 18). When we add the edges of Ai n H back into the graph and 
perform the r-transformation, it is easy to see that the parity of the number of connected 
components is preserved. In the upper case, there are an odd number of path components; 
in the lower case, there are an even number. 

Up to this point, we have shown that C(A) is an invariant of H a t> c . To calculate it, 
we compute the number of edges in the matching corresponding to the empty 3D Young 
diagram. This is easy: there are ab horizontal edges, ac edges of slope v^3, and be edges of 
slope — v3, and so there are ab + bc + ca edges in all. □ 

Lemma 9. If X is a 2-factor on H a ^ c and /1 £ \I/ -1 (A), then u> p ,i,i,i(A) = w p {[i). 



Proof. Comparing the weightings it> p ,i,i,i and w p in Figure 19, it is clear that if e is an edge 
of H a b c and e' is either of the two edges in ?/> _1 (e), then ti> Pl i,i,i (e') = w p (e). The statement 
of the lemma follows immediately. □. 

6. A SPECIALIZATION OF THE Z 2 X Z 2 PARTITION FUNCTION 

As before, w p>q ^ s denotes the Z 2 x Z 2 weighting, and w p = w PtP ^ p denotes the monochro- 
matic weighting. We will be working with various specializations of these weightings, and 
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Figure 18. The two possible connectivities of A^ U A2, and their behaviour 
under r-transformations 





it will be convenient to write, for example, for the specialization p 1— > —p, 

q,r, s 1— > —1. /usr/XHR6/bin/gvim 

Although our tool of choice is weighted matchings on hexagon meshes, our object of study 
is the Z2 x Z2 weighting for 3D Young diagrams. We shall use w p ^ rjS {'K) to denote the 
Z 2 x Z 2 -weight of the 3D diagram n, and w p to denote its monochromatic weight. Recall 
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that if A is the matching corresponding to ir, then 



Wp,q,r,s( n ) 



w P,g,r,s (empty 3D diagram) 
and we can normalize the monochromatic weight in the same fashion. 

Definition 10. The Z2 x Z2 partition function for H a b c is 



Z a > b ' c (p } q } r } s) = J2 



7r 3d diagram 
inside axbxc box 



The monochromatic partition function for H a ^ fi is 

Z"A°(p)= J]) w p (7r) 

tt 3d diagram 
inside axbxc box 

Note that lim^^^ Z a < b > c = hm a)6)C _» 00 Z 2a > 2b > 2c = Z %2 ^ 2 . Thus the following implies Q: 
Theorem 11. 

Z 2a ' 2b ' 2c (p, -1, -1, -1) = (Z a ' b ' c (-p)) 2 . 

Proof. To begin with, let us work with matchings. Observe that for any matching //, we 
have w p -1-1-1(11) = t0- p ,i,i,i(ju) ■ W-i-i-i-i(fjL). Therefore, 



Y w P,-h-i,-M = Y Y w -v,^M w -h--*--i-M 

fi matching A 2-factor fi^^Sf~ 1 (X) 

on H 2a ,2b,2c on H a,b,c 

(5) = Y w -p( x ) Y w -i,-i-i,-M 

A 2-factor /U6*- 1 (A) 
on -ff a ,;,, c 

^ _ ^_^ab+bc+ca ^ 2#{ closcd loops in A}^ 

A 2-factor 
on H a>b>c 

(7) =( _ 1) «N*+«a £ W_ p (\) W _ p ( V ) 

X,n 1-factors 

( V 



-^\ab+bc+ca 



Y w -p(/ i ) 

fj, matching 
\ on H ajb>c J 



In the above sequence of steps, Equation ^ uses Lemma |9j Equation ^ uses Lemma [7J 
Now, let us normalize both sides by dividing by the w p _i weight of the matching 

associated to the (2a, 26, 2c) empty 3D diagram. Let us call this matching M. Call the 

(a, b, c) empty 3D diagram N. 

On the left side, we get Z a,b ' c (p, q, r, s). For the right side, we observe that 

w p ,- ll _i ) -i(M) = iw_i,_i_ a ,_ 1 (A0 ■ u>- P ,i,i,i( M ) 
= (_i)^+- . (w- p (N)) 2 , 
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again using Lemma |9j Therefore, the right-hand side normalizes to 
and we are done. 



□ 



7. Conclusions 

It is obvious to ask whether this method might be used for its intended purpose, namely 
to prove Equation [2} Unfortunately, the answer is no. Our method is clearly agnostic as to 
the large-scale shape of the graph; it works on any suitable subset of the hexagon lattice. 
However, computer calculations show that the analogue of (|2j) on i^2a,2f>,2c does not admit a 
product formula; rather, if one counts 3D Young diagrams fitting inside a 2a x 2b x 2c box, ac- 
cording to the Z 2 x Z 2 weighting, one typically gets (1 —p) x (a large irreducible polynomial); 
the product formula only appears as a, b, c — > oo. 

It should be possible to apply this method to other graphs and other weights. We have not 
investigated this possibility in any great depth, but in order for the method to work without 
serious modification, it seems that the graph ought to have propeller-like vertices, at which 
several hexagonal faces meet; furthermore, one might expect the graph to have some degree 
of self-similarity at different scales. 
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